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Abstract. We consider the Ising model on a Cayley tree of order two with nearest 
neighbor interactions and competing next nearest neighbor interactions restricted to 
spins belonging to the same branch of the tree. This model was studied by Vannimenus 
and found a new modulated phase, in addition to the paramagnetic, ferromagnetic, an- 

tiferromagnetic phases and a (H — I ) periodic phase. Vannimenus's results based on 

the recurrence equations (relating the partition function of an n— generation tree to the 
partition function of its subsystems containing (n — 1) generations) and most results are 
obtained numerically. In this paper we analytically study the recurrence equations and 
obtain some exact results: critical temperatures and curves, number of several phases, 
partition function. 
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1. Introduction 

The model considered by Vannimenus [20] consists of Ising spins (a = ±1) on a Cayley 
tree of branching ratio 2, so that every spin has three nearest-neighbor (NN). Two kinds 
of bonds are present: NN interactions of strength J\ and next-nearest-neighbor (NNN) 
interactions J2, these being restricted to spins belonging to the same branch of the tree. 
The phase diagram described by Vannimenus contains a modulated phase, as found for 
similar models on periodic lattices, but the multicritical Lifshitz point is at zero tempera- 
ture. The variation of the wavevector with temperature in the modulated phase is studied 
in detail, it is shown narrow commensurate steps between incommensurate regions. The 
behavior of the coherence length near the different transitions is also analyzed. 

The Vannimenus's model was then generalized in many directions: 

In [5] it was considered a model with the competing NN and NNN interactions Ising 
model on a Cayley tree but in their case it is allowed for all interbranch NNN interactions 
on the coordination number three which was discussed earlier in [9] and it were obtained 
in addition to the expected paramagnetic, ferromagnetic and antiferromagnetic phases, 
an intermediate range of J2/J1 < values where the local magnetization has chaotic 
oscillatory glass-like behavior. 

Another generalization is due to Mariz et al [10] these authors studied the phase dia- 
gram for the Ising model on a Cayley tree with competing NN interactions J\ and NNN 
interactions J2 and J3 in the presence of an external magnetic field. At vanishing tempera- 
ture, the phase diagram is fully determined, for all values and signs of J2/J1 and J3/J2; in 
particular, it was verified that values of J3/J2 high enough favor the paramagnetic phase. 
At finite temperatures, several interesting features (evolution of reentrances, separation of 
the modulated region into two disconnected pieces, etc.) are exhibited for typical values 
of J2/J1 and J3/J2. 
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The next generalization is considered in [TT], where the lattice spin model with Q- 
component discrete spin variables restricted to having orientations orthogonal to the faces 
of Q-dimensional hypercube is considered on the Cayley tree (Bethe lattice). The par- 
tition function of the model with dipole-dipole and quadrupole-quadrupole interaction is 
presented in terms of double graph expansions. By analyzing the regions of stability of 
different types of fixed points of the system of recurrent relations (which is generalization 
of the Vannimenus's equations), the phase diagrams of the model are plotted. For Q < 2 
the phase diagram of the model is found to have three tricritical points. 

The next generalizations are considered in [2], [3] and These authors have studied 
the phase diagram for Potts model on a Cayley tree with competing NN interactions J\, 
prolonged NNN interactions J p and one level NNN interactions J Q . In |4j the Potts model 
with J Q 7^ is considered and it is shown that for some values of J Q the multicritical Lifshitz 
point be at non-zero temperature and proven that as soon as the same-level interaction 
J is nonzero, the paramagnetic phase found at high temperatures for J Q = disappears, 
while Ising model does not obtain such property. 

But most results of the above mentioned works are obtained numerically. Thus it is 
natural to try to get some these results by an analytical way. 

In this paper we consider the same model which was considered by Vannimenus (not 
its generalization) and study its phases by an analytical argument. Here we shall combine 
analytical arguments of papers [5], [H], [E], |19j . 

2. Preliminaries 

The Cayley tree (Bethe lattice [1]) F k of order k > 1 is an infinite tree, i.e., a graph 
without cycles, such that from each vertex of which issues exactly k + 1 edges. Let 
T k = (V, L), where V is the set of vertices of T k , L is the set of edges of T k . Two vertices 
x and y are called nearest neighbors (NN) if there exists an edge I G L connecting them, 
which is denoted by / = (x,y). A collection of the pairs {x, x\), {xa-i, y) is called a 
path from x to y. Then the distance d(x,y),x,y G V, on the Cayley tree, is the number 
of edges in the shortest path from x to y. 

For a fixed x° G V we set 

n 

W n = {xe V\d{x, x°) = n}, V n =\J W k . 

k=l 

Denote 

S(x) = {y£ W n+1 : d(x, y) = 1}, x £ W n , 
this set is called a set of direct successors of x. 

The vertices x and y are called next-nearest-neighbor (NNN) which is denoted by )x, y{, 
if there exists a vertex z G V such that x, z and y, z are NN. We will consider only prolonged 
NNN )x,y{, for which there exist n such that x G W n and y G Wn+2> this kind of NNN is 
considered in [20] and [3]. 

We consider Ising model with competing NN and prolonged NNN interactions on a 
Cayley tree where the spin takes values in the set $ := {—1, 1}, and is assigned to the 
vertices of the tree |20| . A configuration a on V is then defined as a function x G V >— > 
o~(x) G the set of all configurations is $> v . The (formal) Hamiltonian of the model is 

(2.1) H(a) = -J l Y, °( x )°(y)- J 2 X>(x)<7(y), 

{x,y)eL )x,y( 

where J\, J2 G R are coupling constants and (x, y) stands for NN vertices and }x,y{ stands 
for prolonged NNN. 
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As usual, one can introduce the notions of Gibbs measure (phase) of the Ising model 
with a competing interactions on the Cayley tree [7J, |18| . [13]. 

The standard approach consists in writing down recurrence equation relating the par- 
tition function 

Z n = ^ exp{-/3H(a n )}, 

of n-generation tree to the partition function Z n _\ of its subsystems containing (n — 1) 
generations. In |20j the partition function Z n of the Hamiltonian fj2. If) is given by 



(2.2) 
Here 



u 



(«) , ..(«) 



+ i4 n) +4 n) 



n > 1. 



,( n ) „,( n ) „,( n ) „(«) 



satisfies the following recurrent equation 



(2.3) 



u 



(n+l) 



.(«) 



!-.(»») 
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(n+l) 



1 (rM n) + 6n 2 n) 



+ 6u 



(ft) 



where a = exp(Ji/3), b = exp(J2/3). 

Consider mapping F : u = (m, «2, M3, M4) G i? + — > F{u) = (u^, u' 2 , u' 3 , u'^) G i? + defined 

by 

r 1 {b- l ui + bu 2 f 

v2 



(6 1 n 3 + bu±) 



Then the recurrent equation f)2.3[) can be written ), n > which in the 

theory of dynamical systems is called trajectory of the initial point under action of 
the mapping F. Thus asymptotic behavior of Z n for n — > 00 can be determined by values 
of limu^' i.e., trajectory of under action of F. In this paper we study the trajectory 
(dynamical system) for a given initial point G R 4 . 



3. Dynamics of F 

3.1. Fixed points. In this subsection we are going to define fixed points, i.e., solutions 
to F(u) = u. 

Denote Fix(F) = {u : F{u) = u}. 

We introduce the new variables a = \/a, V{ = y/Hi, i = 1,2,3,4. Then the equation 
F(u) = u becomes as 

v 1 = a {bv\ + 
v 2 = or 1 {bv\ + 
V3 = or 1 {\r l v\ + bv%) 
v 4 = a (b~ 1 vl + bv\) . 



(3.1) 
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Lemma 3.1. If a vector u is a fixed point of F then u G Mi = {u = (u\, 112, U3, u^) G R\ : 
ui = u 4 , u 2 = u 3 } oru G M 2 = {u = (u±, u 2 , U3, U4) G R% : yju% + y/ul = ^{^/U2 + y/uz)} , 
where tp{x) = 

Proof. From (|3.ip we get 

(vi - Vi)[ab(vi + U4) - 1] + (i>2 - V3)[o!& -1 («2 + V3)] = 
(t?i - V4)[(a&) _1 ( v l + Vi)} + (u 2 - i>3)[l + a _1 6(w2 + V3)] = 0. 



(3.2) 



If «i = f4 (resp. «2 = 1*3) from the second equation of (13.2j) we get i>2 = t>3 (resp. i>i = U4). 
Thus v\ = v± if and only if V2 = V3. Assume now v\ / V4 and V2 / W3 then system (|3.2I) 
can be reduced to the equation 

(3.3) (6 2 - 6" 2 )(ui + v A ){v 2 + v 3 ) + ab(vt + u 4 ) - a~ l b{v 2 + v 3 ) - 1 = 0. 

The equation (|3.3p gives i>i + W4 = tp( v 2 + ^3)- E~J 

Let us first study fixed points of F which belong in M%: the condition u\ = U4, U2 = U3 
reduces the equation F(u) = u to the following equation 

1 + 6 2 x N 2 



(3.4) x = f{x) = a 2 

where x = Denote 



b 2 + . 



"2 



d = a - 2 b- 6 , b = b 4 , y = b 2 x. 
The following lemma gives full description of solutions to ([37 

Lemma 3.2. (Cf. Lemma 10.7 in [13] ). Equation j3.^[ ) /ias a unique positive, stable 
solution if b < 9. If b > 9, then there exist v\{b) and ^2(6) suc/z i/mi i/ie conditions 
< ^1(6) < z/2(6) are satisfied and equation {3.$ has three solutions, 
and X3 are stable and x\ is unstable, if v\(h) < a < v^ib) and has two solutions, x\, x\, 
x* is unstable (saddle) and x\ is stable, if a = v\(b) or a = ^2(6)- In this case, we have 

Vi(b) = — 



Vi \b + yi. 

where y\ and 2/2 are the solutions of the equation y 2 + (3 — b)y + 6 = 0. 

Now we shall give some argument to find fixed points of F which belong in M%. Again 
use variables Vi, i = 1, 2, 3, 4, assume V2+V3 = C, with C > max{0, b -z b _ b i }■ Using Lemma 
13.11 we get V3 = C — V2 and V4 = tp(C) — v\. Then from the first equation of (|3.ip we 
obtain V2 = \fb{a~ 1 v\ — bv 2 ) and from the second equation of (13. ip we have ^4(^1) = 
with a polynomial P4 of degree 4, coefficients of which depend on parameters a, 6 and C. 
Thus a quartic equation can be obtained. Such an equations can be solved using known 
formulas (see http://mathworld.wolfram.com/QuarticEquation.html), since we will have 
some complicated formulas for the coefficients and the solutions, we do not present the 
solution here. 

3.2. Periodic points. A point u in R 4 ^ is called periodic point of F if there exists p so 
that F p (u) = u where F p is the pth iterate of F. The smallest positive integer p satisfying 
the above is called the prime period or least period of the point u. Denote by Per p (i ? ) the 
set of periodic points with prime period p. 

Note that the set M\ is invariant with respect to F i.e., F(M%) C M\. In this subsection 
we shall describe some periodic points of F which belong in M\. 
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Let us first describe periodic points with p = 2 on Mi, in this case the equation 
F(F(u)) = u can be reduced to description of 2-periodic points of the function / de- 
fined in (|3.4p i.e., to solution of the equation 

(3.5) f(f(x))=x. 

Note that the fixed points of / are solutions to (|3.5p . to find other solutions we consider 
the equation 

f(f(x))-x _^ 
f(x)-x 

simple calculations show that the last equation is equivalent to the following 

(3.6) 6 4 (1 + a 2 b 2 ) 2 x 2 + a 2 {b 8 + 2(a~ 2 + a 2 )6 6 + 46 4 - l}x + b\a 2 + b 2 ) 2 = 0. 
The equation has two positive solutions iff B < and D > where 

B = a 2 {b 8 + 2(a~ 2 + a 2 )b & + 46 4 - 1}, D = B 2 — {2b\a 2 + 6 2 )(1 + a 2 b 2 )) 2 . 
We have 



' >0, if b>\/y/2-l or b < vV2 - 1, a 2 G (0, b~) U (b+, + 



S = < 



oo 



0, if6<Vv / 2-l, a 2 = b"b+ 



, < 0, if 6 < V\/2- 1, a 2 G (b-,b+) 



where 



(3.7) 



where 



. _ 1 - 46 4 - b 8 ± (1 - 6 4 ) v / (6 4 -l) 2 -46 4 

~ 456 ; 

D = -a 2 (6 4 - l) 2 (46 6 a 4 + (36 s + 66 4 - l)a 2 + 46 6 ) = 

<0, if \/3 _T <6, 6/1 or 6 < y/¥^,a 2 G (0,6~) U (6+,+oo) 
0, if 6 = 1 or 6 < ^/F 1 , a 2 = 6~, 6+ 
[ >0, if b< VF T ,a 2 G (6", 6+) 



1 - 36 8 - 66 4 ± ^(o 4 - 1) 3 (96 4 - 1) 



One can check that V3 -1 < v \71 — 1, and for 6 < a/3" 1 one has b < 6„, and 6+ < b + . 
Thus we have proved the following 

Lemma 3.3. The solutions to 13. 5\) which are different from fixed points of f are vary as 
follows: 

1) 7/v / 3- r < b, b^ 1 or b < V^, a 2 G (0, b~) U (£>+, +oo) i/ien tfie equation [3~M has 
no positive solution. 

2) If b = 1 or b < V3 _1 and a 2 = b~,bf then the equation \3. 6\) has unique positive 



Solution X\— 2fe4(l+a 2 fe 2 ) 2 ' 

3) If b < V3 _1 , a 2 G (0^,0^) i/ien i/iere are too positive solutions x± = 2 b^{i+a, 

(ESP- 



B±^D 



to 



The following lemma gives useful properties of the function /. 
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Lemma 3.4. 1) If b > 1 then the sequence x n = f{x n -\), n = 1, 2, ... converges for any 
initial point xq > 0, where f is defined in {3.$ - 

2) If b < 1 then the sequence y n = g{y n -\), n = 1,2,... converges for any initial point 
y > 0, where g(x) = f{f{x)). 

Proof. 1) For b > 1 we have f'(x) = 2a? (6 4 — 1) /ffi^a > i.e., / is an increasing 
function. Here we consider the case when the function / has three fixed points x*, i = 
1,2,3 (see Lemma 13.21 This proof is more simple for cases when / has one or two fixed 
points). We have that the point x* 2 is repeller i.e., f'(x*£) > 1 and the points x*,x% 
are attractive i.e., fix*) < 1, f'(x^) < 1. Now we shall take arbitrary xq > and 
prove that x n = /(x n _i), n > 1 converges as n — > oo. Consider the following partition 
(0,+oo) = (0,x|)U{^}U(^,4)U{4}U(4,a;|)U{a;|}U(a;|,+oo). For any x G (0,x|) 
we have x < /(x) < a;*, since / is an increasing function, from the last inequalities we get 
x < f(x) < f 2 (x) < f(x\) = x* iterating this argument we obtain f n (x) < f n {x) < x*, 
which for any xq G (0,x*) gives x n _i < x n < x\ i.e., x n converges and its limit is a fixed 
point of /, since / has unique fixed point x\ in (0, x*] we conclude that the limit is x*. 
For x £ (x*,X2) we have x\ > x > f(x) > x*, consequently x n > x n+ \ i.e., x n converges 
and its limit is again x\. Similarly, one can show that if xq > x\ then x n — > x% as n — > oo. 

2) For b < 1 we have / is decreasing and has unique fixed point x\ which is repelling, 
but g is increasing since g'{x) = f'{f{x))f'{x) > 0. By Lemma 13.31 we have that g 
has at most three fixed points (including x\). The point x\ is repelling for g too, since 
g'{x\) = f\f{x\))f'{x\) = (f'(xi)) 2 > 1. But fixed points x_, x + (see Lemma [373]) of g 
are attractive. Hence one can repeat the same argument of the proof of the part 1) for 
the increasing function g and complete the proof. □ 

Lemma 13.31 shows that if b > 1 i.e., J2 > then there is no any 2-periodic trajectory of 
F on Mi. Moreover, the following lemma says that if J2 > then there is no any periodic 
trajectory on Mi. 

Lemma 3.5. 1) If J2 > then for any p > 2 the equation F p (u) = u has no solution 
u £ Mi \Fix(F). 

2) If J2 < then for any p > 3 the equation F p (u) = u has no solution u G Mi \ 
(Fix(F) UPer 2 (F)). 



U ( n +P) = u\ n \ i = 1, 2; n = 0, 1, 2, ... . Since the set Mi is invariant with respect to F, we 



Proof. 1) Assume there is a solution G Mi\Fix(i ? ) then we get p— periodic trajectories 
,.( n +p) . 

obtain 

1 

T\_ 1 — - 

,( n +P) „,(") 



C «+p = ( n+ „l = "7^7 = x n = f U (xo). 

Un 



This is a contradiction, since by Lemma 13.41 we have x n is not periodic. 
2) Assume there is a solution vS® 1 G Mi \ (Fix(i ? ) U Pera(-F)) then we have 

(2n+2p) (2n) 

1 1 f2n { \ 

Vn+p = X 2n +2p = (2n+2 p) = ~~ (2n) ~ Vn ~ J WO)' 

This is a contradiction, since by Lemma 13.41 we have y n is not periodic. □ 



4. Exact values 
idition (with 7 

the recurrence equations (|2.3p and observes their behavior after large number of iterations 



Starting from random initial condition (with 7^ and 7^ ^3°^), one iterates 
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In the simplest situation a fixed point u* = (u* ) G R\ is reached. It corresponds 

to (see [20]): 

a paramagnetic phase if u* G Mi; 

a ferromagnetic phase if u* G M2. 

If the iterations give a cyclic (periodic), say with period p, sequence then the correspond- 
ing phase is called (p-) commensurate phase. Finally, the system may remain aperiodic, 
which corresponds to an incommensurate phase. 

The condition b > 9 of Lemma 13.21 requires that J2 > 0. Denote 

T c = % J 2 >0. 
m3 

Lemma 13.21 gives the following 

Theorem 4.1. If T > T c then the model \2. 1\) has unique paramagnetic phase; If T < 
T c then there are exactly three (resp. two) paramagnetic phases if (J±, J2) is such that 
6 3 v / z7j" < a -1 < b 3 ^2 (resp. aT 1 = b 3 y/v{ or a -1 = b 3 y/i^). 



For the condition b < V3 1 of Lemma 13.31 we need to condition J 2 < 0. In this case we 

ve T c = =ff, J 2 < 0. 

From Lemma 13.31 we get the following 



Theorem 4.2. IfT>T c then the model $2.1\) (on M\) has unique 2- commensurate phase; 
If T < T c then there are exactly two (resp. one) 2- commensurate phases if a 2 G (b~,bf) 
(resp. a 2 = b~ or a 2 = bf ). 

For a fixed temperature T = /3 _1 < T c we have two critical curves a 2 = b^ i.e., on terms 
of J\ and J2 < they are given by the following explicit relations 

J t = _L Ln - 3e 8J ^ _ 6e 4J 2/3 ± ^/( e 4J 2 /3 _ i)3( 9e 4j 2 /3 _ ^jj _ 3 j 2 

Using Lemma 13.31 and formula (I2.2p we can get explicit formulas for the sequence of 
periodic partition functions: 

Z n = Z ri (y) = 2a- 2 ' z x 

"Kb + + + 1) 2 ) ~ § + « § (f (6y + i) 2 + J(lf + ^) 2 )" § ) 2 . - = 2m 
+ to) 2 + a(4 + jfe) 2 )" 1 + J (* (6/(y) + + + b)2 )" f ) 2 ' n = 2m 

where m = 0, 1, 2, . . . ; y is one of xi, x_, x+ defined in Lemma [3.31 and function / is given 
in (I3T41 . 

It is easy to see that if x is a fixed point of / then corresponding fixed point of F 
has the form u*(x) = (u*(x), u^x), u^x), with -uf(x) = a _1 (6 + (6x) -1 ) -2 and 

^2(3:) = a(6 + b~ 1 x)~ 2 . If y is a fixed point of g then corresponding 2-periodic point of F 
has the form u per (y) = {u\ ev {y), u\ ev {y), u\ ev {y), u\ ev {y)) with 

uT\y) = a-V 3 (ab(b + {by)^) 2 + (ab)-\by^ + b- 1 ) 2 )^ 3 , 
Lemma 13.51 gives 
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Theorem 4.3. The model \2. 1\) (on M\) has uncountable set S of incommensurate phases 
fj> u , where u £ M\ \ (Fix(F) U Per 2 (F)). Moreover the set of incommensurate phases can 
be classified to (uncountable) subsets 

S x = {fi u : u G Mi \ (Fix(F) U Per 2 (F)) with lim F n (u) = u*(x)}, 

where x is an attractive fixed point of f and 

SP er = {n u : u € Mi \ (Fix(F) U Per 2 (F)) with lim F 2n (u) = u pcT (y)}, 

where y is an attractive fixed point of g. 



5. Concluding remarks 

Usually, to describe phases (Gibbs measures) of a given Hamiltonian on a Cayley tree 
one has correspondence between Gibbs measures and a collection of vectors (real numbers 
in some particular cases) {h x , x £ V}, which satisfy a non- linear equation (see for example, 
[5], [6], [32]-[l6], [19]). The recurrent equation (12. 3p considered in this paper (which was 
obtained in |20j) describes a vector function {u^ n \n £ N} which is a particular case of 
the above mentioned function h x obtained as h x = vS n > if x € W n i.e., depends only on 
number of the generation set where belongs x but not on x itself. Thus the solutions to 
the recurrent equation (|2.3p do not fully describe phases of the model (|2.ip . But deriving 
of the functional equation for h x corresponding to the Hamiltonian (|2.ip is also difficult, 
since there is prolonged NNN interaction. Such model can be also studied by a contour 
argument (see [T7] and the references therein), but this argument does not give exact 
solutions, in general. 

By a process of iteration, for the model (|2.ip Vannimenus found a new modulated phase, 
in addition to the expected paramagnetic and ferromagnetic (fixed point) phases and a 

(+ H ) periodic (four cycle antiferromagnetic phase, which consisted of commensurate 

(periodic) and incommensurate (aperiodic) regions corresponding to so called "devil's 
staircase". In this paper, using theory of dynamical systems we have analytically proved 
many above mentioned results, i.e., the following exact results are obtained: 

Paramagnetic phase: The exact critical temperature and exact critical curves are found. 
It is proven that the number of the paramagnetic phases can be at most three. (Theorem 

HH). 

Ferromagnetic phase: We reduced description of such phases to a quartic equation (i.e., 
solution of the equations on M 2 ). But we were not able to study the periodic solutions on 
M 2 . 

Commensurate phase: The exact critical temperature (which is obtained from the crit- 
ical temperature of the paramagnetic phase by replacing J 2 with — J 2 ) and exact critical 
curves are found. On the set Mi it is proven that the number of the 2-commensurate 
phases can be at most two and there is not p-commensurate phases if p > 3 (Lemma 13.51 
Theorem 14. 2p . We also described exact values of periodic partition functions. 

Incommensurate phase: We proved that the model has uncountably many such phases. 
Moreover we classified them in two classes: the first class contains the phases which are 
"asymptotically fixed" (set S x ); the second class contains the phases which are "asymp- 
totically periodic" (set Sy C1 ). Note that for the usual Ising model with external field on 
Cayley trees such infinitely many phases are known (see [7j, p. 250). 
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